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AN ANALYSIS AND PROCEDURE FOR DETERMINING SPACE 
ENVIRONMENTAL SINK TEMPERATURES WITH 
SELECTED COMPUTATIONAL RESULTS 


Albert J. Juhasz 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


SUMMARY 

The purpose of this report was to analyze the heat-transfer problem posed by the determination of spacecraft 
temperatures and to incorporate the theoretically derived relationships in the computational code TSCALC. The 
basis for the code was a theoretical analysis of the thermal radiative equilibrium in space, particularly in the Solar 
System. Beginning with the solar luminosity, the code takes into account these key variables: the spacecraft-to-Sun 
distance expressed in astronomical units ( AU), where 1 AU represents the average Sun-to-Earth distance of 
149.6 million km; the angle (arc degrees) at which solar radiation is incident upon a spacecraft surface (ILUMANG); 
the spacecraft surface temperature (a radiator or photovoltaic array) in kelvin; the surface absorptivity-to-emissivity 
ratio aft with respect to the solar radiation and (aJt) 2 with respect to planetary radiation; and the surface view' factor 
to space F. 


INTRODUCTION 

For the thermal design of spacecraft radiators, it is necessary to determine a space background temperature, or 
equivalent space sink temperature, to w hich radiators reject their design heat load. This is especially true for radia- 
tors that operate at the relatively low temperature range of 300 to 500 K, for which early assumptions of 0 to 3 K for 
the sink temperature introduced serious errors in the determination of the required radiator area. This equivalent 
space sink temperature is not the temperature of the vacuum surrounding a spacecraft; rather, it is an equilibrium 
temperature that a passive radiating surface assumes by exchanging thermal radiation energy with the space environ- 
ment without any onboard (internally generated) heat load that has to be rejected to space. This equilibrium tempera- 
ture, a function of the spacecraft-to-Sun distance and the angle at which the surface intercepts solar radiation, is 
greatly influenced by the characteristics (key variables) of the radiating surface. These are the view factor to space F 
and the a/e ratio, which expresses the surface absorptivity of solar radiation to the emissivity at its final equilibrium 
temperature when heat is radiated to space, usually in the infrared region of the spectrum. 

Early reports (Goldman and Singer, 1957) assumed for radiator surfaces gray-body characteristics that gave aft 
values near unity. The work of Rittenhouse and Singletary (1968) showed that special metal oxide coatings such as 
Z-93 could be applied to surfaces to lower the amount of solar energy absorbed while keeping infrared emissivity 
high. With the achievement of a ft values near 0. 1 for these coatings, the space equivalent sink temperatures have 
been effectively lowered, permitting more heat to be rejected by a radiator with a given surface area and an effective 
radiating temperature. 

The purpose of this report was to derive the equilibrium temperatures that actual spacecraft radiating surfaces 
would experience from the combined effects of solar and planetary radiation and to determine the amount of onboard 
heat to be rejected. The derivation is based on a new definition of the term “space sink temperature" and differs from 
earlier definitions, such as that of Gordon (1982), who defined the temperature of space as the equilibrium tempera- 
ture that a small black sphere would experience. 

In contrast to Gordon’s (1982) definition of the space sink temperature, the theoretical derivations in this report 
are based on the equilibrium temperature that a radiating surface having a given aft ratio would achieve in space. 

The conditions are that the surface does not have internal heat to be rejected but is exposed to a given incident radi- 
ant energy flux from the Sun (or star) and planetary surfaces, part of which is reradiated to the space background. 

The amount of thermal energy that reradiates into space is determined by the key variables, such as the aft ratio 
and the view factor to space. Two additional statements can be made regarding the authors’ definition of space sink 
temperature: 
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1. If the incident radiant energy flux is reduced to zero, the equilibrium sink temperature for a surface without 
internal heat generation would drop to 0 K. Note, however, that this condition cannot be realistically achieved, even 
in interstellar space where thermal equilibrium is attained at about 3 K. 

2, For a radiating surface that is rejecting an internal heat load, the average equilibrium surface temperature 
increases to a value consistent with the higher temperature value in the Stefan-Boltzmann equation. Substituting the 
sink temperature, as expressed by the new definition, for the lower temperature value in the same equation results in 
the radiated heat being equal to the internal heat load that must be rejected by the surface. 


SYMBOLS 

A albedo, fraction of incident radiation reflected from a planetary surface 
A RS radiator surface area, m 2 

A s radiating surface sail area, where A RS = A^F is total radiating surface 

d mean distance between the Earth and Sun, AU 

F view factor to space 

/ energy flux (incident energy per second per unit area), W/m 2 

fj individual radiation fluxes, including solar radiation, reflected radiation, and flux due to internal heat 
generation, W/m 2 

fg heat flow’ per unit cross-sectional area, W/m 2 
H orbital altitude, m 

L solar energy generation rate, or luminosity, W 
Q radiative heat transfer, W 

Q w onboard heat generation, W 

R planetary radius (6378 km for Earth), km 

S solar constant 

T e thermal equilibrium at temperature, K 

T ES equilibrium temperature for rotating satellite, K 

T r arbitrary radiating temperature, K 

T ^ equivalent space sink temperature, K 

hot and cold temperatures, respectively, K 

a absorptivity, fraction of incident radiation absorbed by the body (for black bodies, a = £ = 1 ; for gray 
bodies, a = e < 1 ) 

aJt absorptivity -to- emissivity ratio 

(ex /ef> infrared planetary radiation absorptivity-to-emissivity ratio 

£ emissivity of radiating surface, where e = 1 for black body and e < 1 for gray body 
0^ angle at which radiator plate intercepts energy from planet that spacecraft is orbiting, deg 

0^ solar illumination angle, deg 

t| surface fin effectiveness (efficiency) 

p reflectivity, fraction of incident radiation reflected from surface of body 
a Stefan-Boltzmann constant, 5.67xl0 -8 W/m 2 K 4 

T transmissivity, fraction of incident radiation transmitted through body 
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ANALYSIS OF RADIATION HEAT TRANSFER IN SPACE 


This section presents a mathematical analysis of radiation heat transfer in space. The derivation of the constitu- 
tive equations and relationships that describe radiant energy exchange in space has been used as the basis for devel- 
oping the computational code TSCALC. 


Derivation of Equations 

Radiative heat transfer . — The radiative heat transfer Q between two bodies at temperatures T { and 7^ where 
T 1 « T j (negligible reradiation), can be expressed according to the Stefan-Boltzmann Law r as 

e = oefA i (r, 4 -T 2 4 ) (1) 

where 

Q radiated heat flow, W 

c Stefan-Boltzmann constant, 5.67xl0~ 8 W/m 2 K 4 

e emissivity of radiating surface, where 8 = 1 for black body and e <1 for gray body 
F surface area view factor between radiating bodies 
A RS radiating surface sail area, m 2 

A s radiating surface sail area, where radiator surface area A RS = A S F is total radiating surface 
T v T 2 hot and cold temperatures, respectively, K 

Objects in the vacuum of space at arbitrary radiating temperatures T R lose heat by radiation to an environmental 
equilibrium, or space sink temperature 7 5 , that is very near absolute zero (-3 K) in interstellar space. However, as 

shown later in the Results and Discussion section, T s can be several hundred degrees kelvin in the neighborhood of 

the inner planets of the Solar System. Hence, to correctly size the areas for space satellite radiators that operate at 
temperatures between 300 and 400 K, the equivalent space sink temperatures must be determined within an accuracy 
of a few tens of degrees. 

For the space radiation case, equation (1) can be rewritten in terms of T R and T s : 

Q = <kFA s (t$-T s 4 ) (2) 

The equivalent space sink temperatures T§ in the neighborhood of planets in the Solar System was calculated as 
follows. Starting with the solar energy generation rate (or luminosity) L, the energy absorption rate per Earth (or 
planet) unit surface area needs to be balanced by the re-emitted energy rate for thermal equilibrium to be achieved at 
T e or T earth . For this analysis, F is set to unity. If the Sun is represented as a heat source at the center of a sphere 
whose radius is equal to the mean distance between the Earth and Sun d , the figure shows that about 1370 W fall on 
each square meter that is perpendicular to the radiating energy flux. The value 1370 W/nrr is referred to as the solar 
constant at 1 AU (astronomical unit), where d = 1.496x10* * m. The value of the solar constant S can be determined 
by dividing the luminosity L of the Sun by the area of the sphere with a radius equal to 1 AU: 

S = -^y (3a) 

4nd 2 


Substituting L with the value of 3.86x1 0 26 W gives 

5 = 1372.5 W/m 2 forlAU (3b) 

and when rounded to three significant figures becomes 

S = 1370 W/ m 2 (3c) 
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A giant sphere, 1 AU in radius, would 
catch all the Sun's radiative energy. 


Figure 1 . — Sun-to-Earth heat radiation case. 

Note that the result shown in equation (3c) represents the incident energy rate per unit Earth surface area that is per- 
pendicular to the solar energy flux. Before we determine the fraction of this energy flux absorbed by the surface, 
some important relationships must be reviewed in the next section. 

Absorption ^reflection, and transmission of radiation . — When electromagnetic radiation impinges on a body, it 
is partially absorbed, partially reflected, and partially transmitted. The relation between the absorbed, reflected, and 
transmitted energy is 


where 


ot + p + T = 1 


(4) 


a absorptivity, fraction of incident radiation absorbed by the body (for black bodies, a = e . = 1 ; for gray bodies, 
a = e < 1) 

p reflectivity, fraction of incident radiation reflected from surface of body 
I transmissivity, fraction of incident radiation transmitted through body 

For opaque bodies such as the planets, transmissivity x = 0. Equation (4) then can be expressed as 

a + p = l (5a) 


Calculation of Planetary Temperatures in Solar System 

Space radiators in planetary orbits are exposed to solar and planetary radiation; therefore, we must first deter- 
mine the equilibrium temperatures at which planets radiate to spacecraft, even if the spacecraft is in the planetary 
shadow and does not receive any solar radiation. 

Planets also reflect part of the incident solar radiation into space where orbiting spacecraft may receive some of 
it in addition to the direct radiation heating by the Sun. The reflectivity of planets is referred to as the albedo and is 
represented by A in equation (5b). Thus, from equation (5a), the absorptivity can be expressed as 

a = 1 - A (5b) 

To proceed with the calculation of equivalent space sink temperature T s at a given planetary orbit distance from 
the Sun, the energy flux, or incident energy per second per unit area/, is defined as 
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( 6 ) 


x _ Energy rate ( W ) 

/ — O 

Area (rrT ) 

For thermal equilibrium at temperature T E to be achieved. 


/emitted x Radiating area = / abs0 rbed x Absorbing area 


(7) 


where (fig. 1 ) 


./emitted — GET e an( ^ ./absorbed . ,2 

And 


(8a) 


Because the Earth is a rapidly rotating planet, once every 24 hr, its entire surface is considered to be at approxi- 
mately the same temperature. Hence, the total surface energy (or heat) emitted will be 


/emitted x4 ^ 2 < 8b > 

where R is the planetary radius (6378 km for the Earth). 

The total heat intercepted by the Earth from the Sun is the heat flux times the projected Earth area on an imagi- 
nary sphere with a radius of 1 AU. This area is represented as nR 2 . Thus, the total energy absorbed is/ absort(ec j x 7 1/?~. 
Substitute the expressions from equation (7) to obtain 


oeT e xAnR 2 

Solving equation (9) for T E > one obtains 


(1- A)LxnR 2 
And 2 


t e 


(1 -A)L 
\6naed 2 


(9) 


(10) 


Substituting the numerical values forL, d< and a yields, after some simplification. 


T e 


= 279 4 


0 - 4 ) 

Ed 2 


(ID 


where T E is in kelvin and d is in astronomical units (Earth = 1 AU). 

If the approximate value of the Earth’s albedo (A ~ 0.3) and the emissivity (£~ 0.8) are substituted in 
equation (11), we obtain the approximate value of 270 K for T E , which will be used to calculate the radiation from 
Earth to the orbiting spacecraft. 

Recall that equation (11) was derived for rapidly rotating planets such as Earth or Mars. For slowly rotating 
planets like Mercury and the Moon, one must take into account that these bodies receive energy over their projected 
(disk) areas and emit energy, not over their full spherical surface areas but only over the same projected areas 
because the remaining surface area is considered to be too cold to radiate a significant amount of energy back to 
space. For such bodies, the thermal equilibrium is thus established when 


oeT e x nR 2 


(1 - A)L x Tt/? 2 

4nd 2 


(12) 
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and 


= 4 H±zA)L 

V 4nG£d 2 


or 


t e 


= 394 4 


(»-*) 

ed 2 


(13a) 


(13b) 


where, as before, the Sun-Earth distance d is expressed in AU. A comparison of equations (10) and ( 1 3a) shows that 
for slowly rotating planets, the equilibrium temperature is higher by a factor equivalent to the fourth root of the pro- 
jected area (i.e., the ratio of sphere surface-to-disk area), namely the fourth root of 4. 

Apply equation (13a), for d = 1 AU, e = 1 , and A = 0.07, to Earth’s moon to obtain 


T e = 394 387 K 

v r 


(13c) 


which is the maximum temperature at the lunar equator at noon. Note that on the dark side (during the 356-hr lunar 
night), the temperature is -70 K. 

Calculation of Spacecraft Temperatures in Planetary Orbits 

In addition to the electromagnetic radiation from the solar spectrum, spacecraft in planetary orbits receive long- 
wavelength infrared radiation from planets they orbit. 

Near- spherical satellites .— For polyhedral spacecraft geometries that are nearly spherical, equation (10) will 
give a good approximation of the equilibrium temperature due to solar radiation. However, the term that expresses 
the absorptivity-to-emissivity ratio for the planetary surface (1 - A)le is replaced by the absorptivity-to-emissivity 
ratio o tie for the radiating surface material. Thus, for near-spherical shapes, the equilibrium temperature for the rotat- 
ing satellite t es becomes 


t es 



(14) 


where the value of the view factor to space F will normally be close to unity but can be lower if the radiating surface 
is partially occulted by projecting parts of the spacecraft. Lowering the value of F will cause t es to increase. 

Flat-plate radiator surfaces . — For flat-plate radiators, in addition to the ot/e ratio, we need to include the solar 
illumination angle 0 5 , the angle at which the radiator plate intercepts the solar flux. Thus, the equilibrium tempera- 
ture T e due to reflected solar radiation can be expressed as 



( 15 ) 


where I is the transmissivity of the atmosphere. If a satellite is outside a planetary atmosphere, as is normally the 
case, x is set to unity. 


N AS A/TM-200 1 -2 1 0063 


6 




Albedo . — Recall that the albedo A is the fraction of incident radiation reflected from a planetary surface; thus, 
for a surface in Earth orbit at altitude //, we have from equation (8a) 


/reflected 



(16a) 


where ./reflected re ^ ecte ^ radiation energy flux per unit area at orbital altitude H (in meters ), and R £ is the 
radius of the Earth (6.378xl0 6 m). Note that the atmospheric transmissivity (x < 1 ) reduces the reflected radiation 
flux. Of this energy flux, the fraction absorbed by a radiating surface is 


./absorbed /reflected BpCX 


(16b) 


where 0^ is the angle at which the radiator plate intercepts energy from a planet that the spacecraft is orbiting. 
Therefore, due to albedo alone, the equilibrium temperature of a satellite's radiating surface is 


ALxsinG^ — 

!_( r e f 




(16c) 


Earth shine, or direct planetary' radiation . — For direct radiation from Earth or a planet at its previously deter- 
mined equilibrium temperature T earth , the radiation heat flux at or near the planetary surface is 

/radiated = aex ^EARTH ( 

and the equilibrium temperature of a radiating surface at orbital altitude //, due to direct Earth (planet) shine, 
is given by 


Te 



(17b) 


where (a/e A is the infrared planetary radiation absorptivity-to-emissivity ratio of the plate radiator surface and x is 
the transmissivity of the atmosphere. In most cases, (a/e) 2 will be close to unity. 

Internal heat generation . — A spacecraft radiator surface normally rejects heat generated onboard. This rejection 
causes a rise in the equilibrium radiating temperature because of the internal heat flux (the heat flow per unit cross- 
sectional area fq): 


f - 

Q FchArs 


(18) 


where t| is the surface fin effectiveness (efficiency) and Q w is the onboard heat generation in watts. In the absence 
of all other effects, this equilibrium temperature would be represented by 


Te 


Qw 

FtGr\Afts 


(19) 
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Combined effects . — -To determine the combined effects on T £ of the heat-transfer mechanisms discussed in 
the previous sections, the individual radiation fluxes f] are added and set equal to aeT^. When solved for T E , the 
following equation expresses the combined effects of solar radiation, albedo, planetary surface temperature radiation 
in the long-w avelength infrared region, and internal heat generation: 



( 20 ) 


RESULTS AND DISCUSSION 

The equations expressing the thermal equilibrium relationships derived in the section Analysis of Radiation 
Heat Transfer in Space were incorporated in the computer code TSCALC. Several options in this code permit the 
user to evaluate the effects of key input variables, such as solar distance d ( AU), the incident angles for solar and 
planetary radiation 0^ (ILUMANG) and 0^ (ILUM2), the absorptivity-to-emissivity ratio a/e (AE) for solar radia- 
tion, the infrared planetary radiation absorptivity-to-emissivity ratio (oc/e) 2 (AE2), and the view factor to space F 
(FV). Note that FV can have a value up to 2 for flat-plate radiators rejecting heat from both sides. If both sides do 
not have a full 271 steradian view of space, the value will be less than 2. For a sphere that receives thermal radiation 
over its projected area but reradiates over the full surface area, FV can have a maximum value of 4. 

An example showing the increase in equilibrium temperatures for a spacecraft surface approaching the Sun 
from the heliopause at 220 AU is shown in table I. Note that the corona temperature of nearly 2000 K at the 
0.02-AU position represents the expected value at the closest approach during the future perihelion mission. A 
carbon-carbon heat shield having an a/e ( AE) of 0.6 and intercepting solar radiation at 25° (ILUMANG) may be 
used to shield the spacecraft. 

Figure 2 graphically depicts the heat flux and temperature information. Note that the temperatures shown are for 
planetary orbit distances, not for the planets themselves. However, planetary temperatures could be computed by the 
code if the proper input values (ILUMANG, FV, and AE) were used. 

Table II shows the effect at 1 AU of the insolation angle, w hich is defined as the angle at w hich a reradiating 
surface without internal heat generation ( Q w = 0) intercepts the solar energy flux. Note that for a solar incident angle 
0 5 of 0°, the code gives a value of 3 K because the solar heat flux has been eliminated and the equilibrium tempera- 
ture is that of interstellar space. 

The information in table III can be used to consider the example of equilibrium temperatures determined prima- 
rily by planetary radiation in the infrared domain for a space structure in geostationary orbit at a 35 876-km altitude. 
By setting the incident angle for solar radiation (ILUMANG) equal to 0.1°, we assume that 99.9 percent of the solar 
heat flux can be eliminated by shielding and insulation. Also, by assuming an Earth-gravity-gradient-stabilized 
structure, one can see how the space sink temperature varies as the direct radiation from the Earth is incident upon a 
surface from 1° to 10° of arc. From the table, we can see that temperatures lower than 77 K (liquid nitrogen boiling 
point at atmospheric pressure) can be maintained if the angle of incident Earth radiation can be held to less than 5°. 


TABLE 1 .—CONDITIONS FOR SPACECRAFT APRO ACHING SUN TO AU = 0.02 


ILUMANG 

(DEG) 

FV 

EPS 

AE 

AU 

Q/A1W/M2) 
(90 DEG.) 

TS(K) 

ORBIT 

25.00 

1.0 

.88 

.60 

220.000 

.03 

18.9 

HLPAUSE 

25.00 

1.0 

.88 

.60 

39.438 

.88 

44.6 

PLUTO 

25.00 

1.0 

.88 

.60 

30.058 

1.52 

51.1 

NEPTUNE 

25.00 

1.0 

.88 

.60 

19.182 

3.73 

63.9 

URANUS 

25.00 

1.0 

.88 

.60 

9.539 

15.08 

90.6 

SATURN 

25.00 

1.0 

.88 

.60 

5.203 

50.70 

122.7 

JUPITER 

25.00 

1.0 

.88 

.60 

3.000 

152.50 

161.6 

ASTRDS 

25.00 

1.0 

.88 

.60 

1.524 

591.18 

226.8 

MARS 

25.00 

1.0 

.88 

.60 

1.000 

1372.51 

279.9 

EARTH 

25.00 

1.0 

.88 

.60 

.723 

2623.26 

329.1 

VENUS 

25.00 

1.0 

.88 

.60 

.387 

9164.15 

450.0 

MERCURY 

25.00 

1.0 

.88 

.60 

.020 

3431265.02 

1979.3 

CORONA 
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Figure 2. — Solar heat flux and spacecraft surface temperatures at various planetary orbit 
distances. 


TABLE II.— SINK TEMPERATURE (TS) AS FUNCTION OF INSOLATION ANGLE AT 1 AU 


ILUMANG 

(DEG) 

FV 

EPS 

AE 

AU 

TS(K) 

QW-WATTS 

.00 

1.5 

.88 

.18 

1.000 

3.0 

.0 

1.00 

1.5 

.88 

.18 

1.000 

84.4 

.0 

2.00 

1.5 

.88 

.18 

1.000 

100.3 

.0 

3.00 

1.5 

.88 

.18 

1.000 

111.0 

.0 

4.00 

1.5 

.88 

.18 

1.000 

119.3 

.0 

5.00 

1.5 

.88 

.18 

1.000 

126.1 

.0 

10.00 

1.5 

.88 

.18 

1.000 

149.9 

.0 

15.00 

1.5 

.88 

.18 

1.000 

165.6 

.0 

20.00 

1.5 

.88 

.18 

1.000 

177.5 

.0 

25.00 

1.5 

.88 

.18 

1.000 

187.2 

.0 

30.00 

1.5 

.88 

.18 

1.000 

195.2 

.0 

35.00 

1.5 

.88 

.18 

1.000 

202.0 

.0 

40.00 

1.5 

.88 

.18 

1.000 

207.9 

.0 

45.00 

1.5 

.88 

.18 

1.000 

212.9 

.0 

50.00 

1.5 

.88 

.18 

1.000 

217.2 

.0 

55.00 

1.5 

.88 

.18 

1.000 

220.9 

.0 

60.00 

1.5 

.88 

.18 

1.000 

224.0 

.0 

90.00 

1.5 

.88 

.18 

1.000 

232.2 

.0 
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TABLE ill.— GEOSTATIONARY ORBIT EQUILIBRIUM TEMPERATURE AT NEAR ECLIPSE 


ILUMANG (DEG) 
Sun Earth 

ALTITUDE 

(K-KM) 

FV 

TR-K 

TS-K 

AE 

AE-2 

AU 

QW-WATTS 

.1 1.0 

35.88 

1.0 

300.0 

59.5 

.18 

.80 

1.000 

.0 

.1 2.0 

35.88 

1.0 

300.0 

64.6 

.18 

.80 

1.000 

.0 

.1 3.0 

35.88 

1.0 

300.0 

68.7 

.18 

.80 

1.000 

.0 

.1 4.0 

35.88 

1.0 

300.0 

72.2 

.18 

.80 

1.000 

.0 

.1 5.0 

35.88 

1.0 

300.0 

75.2 

.18 

.80 

1.000 

.0 

.1 6.0 

35.88 

1.0 

300.0 

77.9 

.18 

.80 

1.000 

.0 

.1 7.0 

35.88 

1.0 

300.0 

80.4 

.18 

.80 

1.000 

.0 

.1 8.0 

35.88 

1.0 

300.0 

82.6 

.18 

.80 

1.000 

.0 

.1 9.0 

35.88 

1.0 

300.0 

84.7 

.18 

.80 

1.000 

.0 

.1 10.0 

35.88 

1.0 

300.0 

86.6 

.18 

.80 

1.000 

.0 


Results of the type illustrated above have been used in the design of radiators for deep space probes (Juhasz, 
et al„ 1999). The computational code can also be used to establish environmental conditions for structures that gen- 
erate solar power in geostationary orbits. 


CONCLUDING REMARKS 

A detailed derivation of the constitutive equations and relationships that describe radiant energy exchange in 
space has been completed and used as the basis for developing a computational code. The equations describe heat- 
transfer mechanisms related to key variables such as solar or stellar radiation, reflected planetary radiation (albedo), 
direct planetary radiation in the infrared portion of the electromagnetic spectrum, and onboard generated heat that 
must be rejected to the space environment. 

To rapidly evaluate a combination of these heat-transfer mechanisms, the comprehensive multioptional com- 
puter code TSCALC was developed and verified against some independent determinations and measurements found 
in the literature. 

Outputs from the code have also been used to determine equilibrium space sink temperatures for deep-space 
probes. 
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